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A three-dimensional lattice of the Josephson junctions with a finite self-conductance is employed 
to model the ceramic superconductors. The nonlinear ac susceptibility and the compensation effect 
are studied by Monte Carlo simulations in this model. The compensation effect is shown to be due 
to the existence of the chiral glass phase. We demonstrate, in agreement with experiments, that 
this effect may be present in the ceramic superconductors which show the paramagnetic Meissner 
effect. 
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One of the most fascinating discoveries in con- 
densed matter physics is the paramagnetic Meiss- 
ner effect (PME) in certain ceramic superconduc- 
tors 0,|2| . The nature of the unusual paramagnetic 
behaviour may be related to the appearance of spon- 
taneous suppercurrents (or of orbital moments) |}) . 
The latter appear due to the existence of 7r-junctions 
characterized by the negative Josephson couplings 
Hfi}. Furthermore, Sigrist and Rice argued that the 
PME in the high-T c superconductors is consistent 
with the d-wave superconductivity Q . This effect is 
succesfully reproduced in a single loop model H as 
well as in a model of interacting junction- loops 

The mechanism of the PME based on the d-wave 
symmetry of the order parameter remains ambigu- 
ous because it is not clear why this effect could not 
be observed in many ceramic materials. More im- 
portantly, the paramagnetic response has been seen 
even in the conventional Nb fj|-|lfH and Al |fl| super- 
conductors. In order to explain the PME in terms of 
conventional superconductivity one can employ the 
idea of the flux compression inside of a sample. Such 
phenomenon becomes possible in the presence of the 
inhomogeneities []l2| or of the sample boundary [[[3) . 
Thus the intrinsic mechanism leading to the PME is 
still under debate [|l]Jl|] . 

Recently Heinzel et al. have shown that the 
PME may be analyzed by the compensation tech- 
nique based on the measurement of the second har- 
monics of the magnetic ac susceptibility. Their key 
observation is that the so called compensation effect 
(CE) appears only in the samples which show the 
PME but not in those which do not. Overall, this ef- 
fect may be detected in the following way. The sam- 
ple is cooled in the external dc field down to a low 
temperature and then the field is switched off. At 
the fixed low T the second harmonics are monitored 



by applying the dc and ac fields to the sample. Due 
to the presence of non-zero spontaneous orbital mo- 
ments the remanent magnetization or, equivalcntly, 
the internal field appears in the cooling process. If 
the direction of the external dc field is identical to 
that during the field cooled (FC) procedure, the in- 
duced shielding currents will reduce the remanence. 
Consequently, the absolute value of the second har- 
monics |%2 1 decreases until the signal of the second 
harmonics is minimized at a field Hd c — H com . Thus 
the CE is a phenomenon in which the external and 
internal fields are compensated and the second har- 
monics become zero. 

The goal of this paper is to explain the CE theo- 
retically by Monte Carlo simulations. Our starting 
point is based on the possible existence of the chiral 
glass phase in which the remanence necessary 
for observing the CE should occur in the cooling 
procedure. Such remanence phenomenon is similar 
to what happens in spin glass. Furthermore, the 
PME related to the CE can also be observed in the 
chiral glass phase j7 16 IT], 
perimental results 



There are several ex- 
which appear to corroborate 
the existence of such a novel glassy phase in ceramic 
high-T c superconductors. 

In the chiral glass phase the frustration due to 
existence of 0- and 7r-junctions (0-junctions corre- 
spond to positive Josephson contact energies) leads 
to non-zero supercurrents [ fl6| . The internal field (or 
the remanent magnetization) induced by the super- 
currents in the cooling process from high tempera- 
tures to the chiral glass phase may compensate the 
external dc field. 

We model ceramic superconductors by the three- 
dimensional XY model of the Josephson network 
with finite self-inductance. We show that in the 
FC regime the CE appears in the samples which 
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show the PME but not in those containing only In- 
junctions. In the zero field cooled (ZFC) regime de- 
creasing the external dc field also gives rise to the 
CE in the frustrated ceramics. Both of these find- 
ings agree with the experimental data of Heinzel et 
a/@. 

We neglect the charging effects of the grain and 
consider the following Hamiltonian |(|f?J 



H = - Jij cos(0i - 9j - A i:j ) + 
<ij> 

P 

p- J2 A * » *v = 2 -r f A ^ » (!) 



<i]> 



where Oi is the phase of the condensate of the grain 
at the i-th site of a simple cubic lattice, A is the fluc- 
tuating gauge potential at each link of the lattice, 
denotes the flux quantum, Jij denotes the Josephson 
coupling between the i-th and j-th grains, £ is the 
self-inductance of a loop (an elementary plaquette) , 
while the mutual inductance between different loops 
is neglected. The first sum is taken over all nearest- 
neighbor pairs and the second sum is taken over all 
elementary plaquettes on the lattice. Fluctuating 
variables to be summed over are the phase variables, 
9i, at each site and the gauge variables, Aij, at each 
link. <j> p is the total magnetic flux threading through 
the p-th plaquette, whereas ^ xt is the flux due to an 
external magnetic field applied along the z-direction, 



HS if p is on the < xy > plane 
otherwise , 



(2) 



where S denotes the area of an elementary plaquette. 
The external field H includes the dc and ac parts and 
it is given by 



H = H dc + H ac cos{u>t) 



(3) 



It should be noted that the dc field is necessary to 
generate even harmonics. 

In the present paper, we consider two models with 
two types of bond distributions. Model I: the sign of 
the Josephson couplings could be either positive (0- 
junction) or negative (7r-junction) and the spin glass 
type bimodal (± J) distribution of Jij is taken. The 
coexistence of 0- and 7r-j unctions gives rise to frus- 
tration even in zero external field and the chiral glass 
phase may occur at low temperatures |l6[] . Model 
II: the interactions Jij are assumed to be 'ferromag- 
netic' and distributed uniformly between and 2J. 



Obviously, there is no frustration in zero external 
field in this model. It has been also demonstrated 
that the PME is present in model I but not in model 

n 0. 

The ac linear susceptibilty of models I and II has 
been studied M by Monte Carlo simulations. It was 
found that, due to the frustration, model I exhibits 
much stronger dissipation than model II in the low 
frequency regime. Here we go beyond our previ- 
ous calculations of the linear ac susceptibility j?]]. 
We study the dependence of the second harmonics 
as a function of the dc field. In this way, we can 
make a direct comparison with the CE observed in 
the experiments 15 1. The second harmonics of a 



similar Josephson network model with a finite self- 
inductance were considered by Wolf and Majhofer 
[p^| . However, these authors dealt with the two- 
dimensional version of model II and the CE has not 
been studied. In this paper we are mainly interested 
in the CE in the frustrated three-dimensional system 
described by model I. 
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FIG. 1. The dependence of \\2\ on hdc for model I 
and model II at T — 0.1. The values of C are chosen to 
be equal to 1, 3 and 10 as shown next to the curves. The 
results are averaged over 20 samples. 

The dimcnsionless magnetization along the z-axis 
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mormalized per plaquette, in, is given by 
1 



Carlo may actually give reasonable results for the 
CE. 



N p (j) 



(4) 



p<xy> 



where the sum is taken over all N p plaquettes on the 
< xy > plane of the lattice. The real and imaginary 
parts of the ac second order susceptibility x^O^O an d 
x'ii 1 ^) are calculated as 



XaM = 

X2M = 



771(f) cos(2a;f) d(u>t) , 
m(f) sm(2ut)d(cjt) , (5) 



where f denotes the Monte Carlo time. The dimen- 
sionless ac field h ac , dc field h^c and inductance C 
are defined as follows 



2vH ac S 



2nH dc S 



L = (2tt/cI >q ) 2 J£. (6) 

The dependence of C on the parameters of the sys- 
tem such as the critical current and the typical size 
of the grains is discussed in [f 17|- 

Our results have been obtained by employing 
Monte Carlo simulations based on the standard 
Metropolis updating technique. While Monte Carlo 
simulations involve no real dynamics, one can still 
expect that they give useful information on the long- 
time behavior of the system. In fact, the amplitude 
of the ac field we use is much smaller than the typ- 
ical energy of the dc part. On the other hand, the 
characteristic time for the sintered samples, which 
are believed to be captured by our model, is of order 
10~ 12 s This time has the same order of mag- 

nitude as a single Monte Carlo step. So the period 
of oscillations chosen in the present work is much 
longer than the characteristic time (see below). For 
such a weak and slowly changing ac field the system 
can be regarded as being in quasi-equilibrium and 
the Monte Carlo updating may be applied. A priori, 
the validity of this approximation is not clear but it 
may be justified by comparing our results with those 
obtained by other approaches to the dynamics such 
as considered in ref. p9[ . For the first harmonics, 
our method and the method of ref. |l9| yield results 
that agree qualitatively. Furthermore, our results 
presented in Fig. 1 for the second harmonics are also 
in a qualitative agreement with the corresponding 
results obtained by solving the equations of motion 
p9|. So one can expect that the standard Monte 
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FIG. 2. The second harmonics of model I obtained 
after field cooling in a dc field hd c — 1 from T = 0.7 to 
T — 0.1. The temperature is reduced in steps of 0.05. At 
the lowest T = 0.1 the dc field used in cooling is switched 
off and the second harmonics are generated by applying 
the combined field (3). The dc field is stepwise reduced 
from hdC = 1 to h^c = 0. The inductance is chosen 
to be equal to £ = 4. The arrows indicate the sense 
of the changes in the dc field. The results are averaged 
over 40 samples and are qualitatively the same as those 
presented in Fig. 1 of Ref. [15]. 

We choose the gauge where the bond variables Aij 
along the z-direction are fixed to be zero. The lattice 
studied are simple cubic with LxLxL sites and free 
boundary conditions are adopted. In all calculations 
presented below, we take L — 8 and u> = 0.001. The 
sample average is taken over 20-40 independent bond 
realizations. X2{w) na s been estimated following the 
procedure in fi^Cfl. Namely, at the beginnig of a 
given Monte Carlo run, we first switch on the field 
(3). Then, after waiting for initial to Monte Carlo 
steps per spin (MCS), we start to monitor the time 
variation of the magnetization, to is being chosen 
so that all transient phenomena can be considered 
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extinct. We set to to be 2 x 10 4 MCS. After passing 
the point t = to, fh{t) is averaged over typically 200 
periods, each period contains ty MCS (tx = 2n/u>). 
The real and imaginary parts of the second order ac 
susceptibility are then extracted via Eq. (5) . We set 
h ac = 0.1, corresponding to w 0.016 flux quantum 
per plaquette. Smaller value of h ac turned out to 
leave the results almost unchanged. 




^dc 

FIG. 3. The procedure to generate the second har- 
monics is the same as in Fig. 2 but for model II. The 
system is cooled in a dc field hdc = 1 from T = 1.4 to 
T = 0.1. The results are averaged over 20 samples. 



The dependence of 1x2 1 = \J{x' 2 ) 2 
on hdc at T = 0.1J is presented in Fig.l. For small 
values of £, the oscillation of \x2\ shows up. Such 
oscillation has been found for the two-dimensional 
superconductors in Ref. and its nature is related 
to the lattice periodicity. Our new observation is 
that the oscillatory behavior is still present in the 
superconductors with 0- and 7r-junctions (model I) 
but to less extent compared to model II. It is clear 
from Fig. 1 that | X2 1 does not decrease at large 
hdc but gets saturated. This is an artifact of the 
assumption that the Josephson contact energies 
are field-independent. The field dependence of J.^ 
should remove the saturation of \x%\ at strong dc 



fields glfy). 

In order to study the difference between model I 
and model II through the CE we have to consider the 
weak field region where the PME may be observed. 
For model I the PME appears clearly for hdc < 1 
0. So the largest hdc we take is 1. In this weak field 
regime there is no periodicity of \xi \ versus hd c which 
may complicate the study of the CE. The chiral glass 
phase is found to exist below a critical value of the 
inductance £ c where 5 < £ c < 7 |l6"| . One has 
to choose, therefore, an £ which is smaller than its 
critical value and in what follows we take £ = 4. 

In this paper we focus on the system size L = 8, 
£ = 4, u = 0.001, and T = 0.1. Our preliminary 
studies show that the qualitative results do not de- 
pend on the choise of the parameters of the system. 

Fig. 2 shows the dependence of second harmon- 
ics on hdc in the FC regime for the superconduc- 
tors described by model I. Our calculations follow 
exactly the experimental procedure of Heinzel et al 
[jl5[ . First the system is cooled in the dc field hdc = 1 
from T = 0.7 down to T = 0.1 which is below the 
paramagnet-chiral glass transition temperature T c ks 

0. 17 |16| . The temperature step is chosen to be equal 
to 0.05. At each temperature, the system is evolved 
through 2xl0 4 Monte Carlo steps. When the lowest 
temperature is reached the dc field used in cooling is 
switched off and we apply the combined field given 
by Eq. (3). We monitor the second harmonics re- 
ducing the dc field from hdc = 1 to zero stepwise 
by an amount of Ahd c — 0.05. \xu\ reaches mini- 
mum at the compensation field h CO m — 0.7 ± 0.05. 
At this point, similar to the experimental findings 
|pl| , the intersection of x'2 and x'2 1S observed. This 
fact indicates that at H CO m the system is really in 
the compensated state. Furthermore, in accord with 
the experiments, at the compensation point the real 
and imaginary parts should change their sign Jl5"| |. 
Our results show that x'2 changes its sign roughly 
at hdc — h CO m- A similar behavior is also displayed 
by x'2 but it is harder to observe due to a smaller 
amplitude of x'i- 

Fig. 3 shows the dependence of the second har- 
monics on hdc in FC regime for model II. The calcu- 
lations are carried out in the same way as for model 

1. A difference is that we start to cool the system 
from T = 1.4 which is above the superconducting 
transition point T s w 0.9 (T s is estimated from the 
maximum of the specific heat (7) for £ = 4 and the 
results are not shown here). The temperature step 
is set equal to 0.1. Obviously, |x2 1 decreases with 
decreasing hdc monotonically. Thus, there is no CE 
because the remanent magnetization does not ap- 
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pear in the cooling process. This result is again in 
accord with the experimental data |lq1 . 
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FIG. 4. The dependence of X2I on hdc obtained in 
the ZFC regime for model I. The solid circles correspond 
to the case when the dc field is decreased from hdc — 1 to 
-0.5. The open hexagons and squares correspond to the 
increase of hdc from zero to 1 and to its decrease from 
zero to -0.5, respectively. The inductance is chosen to be 
equal to L = 4. The sense of changes of the dc field is 
marked by the arrows. The results are averaged over 25 
samples. The compensation field in the case when the 
field is decreased is h CO m = 0.15 ± 0.05. 

We now turn to the ZFC regime. The experiments 
[PH show that no CE can be expected if after the 
ZFC procedure one increases the dc held. However, 
if the held is decreased a remanent magnetization 
is developed and the CE appears [|l5|. The results 
of our simulations for the ceramic superconductors 
described by model I are shown in Fig. 4. As in 
the FC regime the system is cooled from T = 0.7 
to T = 0.1 but without the external held. Then at 
T = 0.1 we apply the held given by Eq. (3) and 
study three cases. In one of them hdc is decreased 
from hd c = 1 to -0.5. The values of \x%\ are repre- 
sented by solid circles in Fig. 4. The CE is clearly 
seen at h com = 0.15 ± 0.05. At this point the real 
and imaginary parts of the second harmonics also 
intersect (the results are not shown). It is not sur- 
prising that the h com in the ZFC regime appears to 
be smaller than in the FC regime. Fig. 4 shows 
also the dependence of \x%\ on the dc held when it 
changes from hdc = to 1 (open hexagons) and from 
hdc = to -0.5 (open squares). Obviously, no CE is 



observed in this case. The results presented in Fig. 
4 qualitatively agree with those shown in Fig. 2 of 
Ref. |5|. 

In conclusion we have shown that the CE may be 
explained, at least qualitatively, by using the chiral 
glass picture of the ceramic superconductors. The 
CE is shown to appear in the chiral glass phase in 
which the PME is present but not in the samples 
without the PME. 
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